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Free-standing Bi 2 Se 3 single crystal flakes of variable thickness are grown using a catalyst-free 
vapor-solid synthesis and are subsequently transferred onto a clean Si ++ /Si0 2 substrate where the 
flakes are contacted in Hall bar geometry. Low temperature magneto-resistance measurements are 
presented which show a linear magneto-resistance for high magnetic fields and weak anti-localization 
(WAL) at low fields. Despite an overall strong charge carrier tunability for thinner devices, we find 
that electron transport is dominated by bulk contributions for all devices. Phase coherence lengths l$ 
as extracted from WAL measurements increase linearly with increasing electron density exceeding 
1 pm at 1.7 K. While l$ is in qualitative agreement with electron electron interaction-induced 
dephasing, we find that spin flip scattering processes limit l$ at low temperatures. 


I. INTRODUCTION 

Topological insulators (TIs) are a new class of mate¬ 
rials m, consisting of an insulating bulk and topologi¬ 
cally protected conducting surface states. These surface 
states are spin polarised and robust against scattering 
from non-magnetic impurities making them interesting 
candidates for future spintronics and quantum comput¬ 
ing devices nu as well as potential hosts for Majo- 
rana fermions [SHU]. Binary Bi-chalcogenides (Bi 2 Se 3 , 
Bi 2 Te 3 ) belong to the class of three-dimensional (3D) 
strong topological insulators with a single Dirac cone at 
the surface m which is experimentally observable by 
angle-resolved photo emission spectroscopy El- In par¬ 
ticular, Bi 2 Se 3 with a single Dirac cone centered in a bulk 
band gap of E g ~ 350 meV is a promising material for 
probing surface states by electronic transport. However, 
the measurement of pure surface states is challenging. 
So far Bi 2 Se 3 crystals are unintentionally n-type doped 
most likely by Se vacancies mm leading to bulk con¬ 
ductivity dominating electronic transport. To increase 
the surface-to-bulk ratio, ultra-thin flakes with thickness 
of the order of 10 nm have been investigated [IS]. 

Thin Bi 2 Se 3 crystals can be produced by mechanical 
exfoliation of bulk material as it is common practice for 
graphene fabrication mm- However, it is much more 
promising to grow thin Bi 2 Se 3 films in-situ which has 
been successfully achieved with molecular beam epitaxy 
(MBE) [2TH23] or vapor-solid synthesis (VSS) in a tube 
furnace [H|2l|25]. 

In this work, we show a catalyst free growth of large 
freestanding Bi 2 Se 3 flakes with a VSS method. Our free¬ 
standing growth approach ensures the synthesis of strain 
free, few-layer single crystal flakes with lateral dimen¬ 
sions up to 25 /im and thicknesses in the range of 6 
to 30 nm, ranking our flakes among the largest Bi 2 Se 3 
single crystals flakes. The high structural and surface 
quality of the Bi 2 Se 3 crystals is verified by Raman and 
by scanning force microscopy. The free standing sin¬ 
gle crystals are ideal for transport studies. We utilize 
a wet chemistry-free process which allows transferring 
these single crystals onto any desired substrate with¬ 


out introducing additional contamination. We studied 
low temperature magneto-transport on a series of Bi 2 Se 3 
crystals of different thicknesses which were transferred 
on Si0 2 /Si ++ substrates. We observe linear magneto¬ 
resistance (LMR) at high H-fields as well as weak anti¬ 
localization (WAL) which both indicate the dominance 
of bulk transport contributions. Electron phase coher¬ 
ence lengths are in the micrometer range, slightly larger 
compared to earlier studies on Bi 2 Se 3 crystals grown di¬ 
rectly on Si0 2 m or by other growth methods [271 EH]. 
We show that electron spin-flip processes limit the phase 
coherence length at low temperatures. 


II. CRYSTAL GROWTH AND 
CHARACTERISATION 

With the goal of gaining high quality thin Bi 2 Se 3 
crystals we applied a catalyst-free vapor-solid synthesis 
method. Most commonly, MBE j29ll3T] is used to grow 
thin films since it offers the growth of extended films of 
rich chemical compositions with excellent thickness con¬ 
trol. Yet it suffers (i) from poly-crystallinity of the film 
inducing strain at the grain boundaries and (ii) from a 
limited number of usable substrates. In contrast to MBE, 
VSS allows the growth of single crystalline platelets on 
a variety of different substrates [32H34j . However, strain 
can still be induced by the growth substrate. Moreover, 
growth catalyst may induce unwanted dopants into the 
crystal. In this work we therefore use a catalyst-free VSS 
method where flakes and ribbons grow free standing on 
the substrate. This offers an interesting pathway for the 
fabrication of high quality devices. Free standing flakes 
are neither strained nor contaminated by the substrate 
material and can be easily transferred onto on a wide 
range of different substrates. 

A standard three zone tube-furnace (Fig. [IJa)) with 
electric heating coils is used for the VSS growth. As a 
source material we place Bi 2 Se 3 crystals [35] in the first 
zone. The Si/Si0 2 growth substrates are placed down¬ 
stream in the second zone. Their exact position was op¬ 
timized through several growth cycles. Prior to growth, 
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the quartz tube was evacuated to 2 mbar with subse¬ 
quent argon flushing for 5 min with 500 seem flow rate 
which is regulated by a digital mass flow controller. After 
cleaning, the growth zone (second zone) is heated up to 
325°C with a constant argon flow of 100 seem to carry 
away vaporized particles. The second zone is kept at 
325°C and 25 mbar for 2 h as it is crucial for the temper¬ 
ature and pressure to be stabilized during growth. Fi¬ 
nally, the actual growth process is executed by heating 
the first zone to 700° C. The source material gradually 
vaporizes and gets carried downstream by an 60 seem 
argon flow. Temperature and pressure were optimized 
to grow large, thin, free standing Bi 2 Se 3 flakes and rib¬ 
bons as shown by scanning electron microscope (SEM) 
images in Figs, [ljb)-(d) and by atomic force microscope 
(AFM) images in Figs. me) -(g). Straight edges with only 
60° and 120° corners indicate single crystalline growth. 
According to AFM measurements the flake thicknesses 
range between 6 nm and 40 nm and lateral dimensions 
can reach up to 25 fi m. AFM images also reveal the 
flake’s surfaces to be stepless confirming a very homoge¬ 
neous layer by layer growth. 

Raman spectroscopy has emerged as an excellent tool 
to probe crystal stoichiometry of Bi 2 Se 3 p6ff38] . The 
Raman spectra of our Bi 2 Se 3 flakes is obtained using con- 
focal Raman spectroscopy with a laser spot diameter of 
around 500 nm at a wavelength of 532 nm. The laser 
spot is precisely positioned on the flakes using a piezo 
stage. In Fig. 1(h) all four characteristic Raman peaks 
of Bi 2 Se 3 are clearly seen at 37 cm -1 , 71 cm -1 , 131 cm -1 
and 175 cm -1 , which correspond to the E g , A\ g , E g and 
A\ g vibrational modes, respectively. The peak positions 
are very close to previously measured Raman peaks of 
stoichiometric Bi 2 Se 3 crystals [36, 39, 40 indicating the 
high crystal quality of our flakes. 

For transport studies, the free standing Bi 2 Se 3 flakes 
are dry-transferred by gently dabbing a clean room cloth 
onto the grown chips and subsequently onto a clean 
Si0 2 /Si ++ substrate. This method does not involve sol¬ 
vents or other liquids which could effect the surface qual¬ 
ity. These substrates are pre-patterned with gold mark¬ 
ers to relocate individual flakes and enable consecutive 
electron beam lithography (EBL). 

The flakes to be contacted are first chosen by optical 
microscopy and further characterized by AFM, which is 
also used to determine exact dimensions. The contacts 
are defined using standard EBL techniques and 5/50 nm 
Cr/Au ohmic contacts. Directly before metal evapora¬ 
tion the contact areas are etched for 15 s by oxygen 
plasma to remove any oxide layer from the Bi 2 Se 3 sur¬ 
face. This step is crucial for low contact resistances. The 
contact geometry of flakes with a high length/width ra¬ 
tio, resembles a fairly good approximation of a Hall bar 
geometry (c/. Figs. j2ja),(b)). With this method no ad¬ 
ditional patterning step is needed, allowing us to keep 
unetched flake edges, as grown in the VSS-process. The 
dimensions of the four investigated devices are summa¬ 
rized in Table [H 
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Figure 1. (color online) (a) Schematic illustration of the 
three-zone oven used for Bi 2 Se 3 sample growth, (b), (c) 
and (d) Scanning electron microscope images of typical ’free 
standing’ Bi 2 Se 3 ribbons and flakes, (e), (f) and (g) AFM im¬ 
ages of grown flakes transferred onto Si0 2 /Si substrate, (h) 
Characteristic Raman spectra of grown Bi 2 Se 3 flakes with dif¬ 
ferent thicknesses. 


III. RESULTS AND DISCUSSION 


Transport measurements were performed in a 4 He- 
cryostat at a base temperature of T = 1.7 K using 
low-frequency lock-in techniques. A superconducting 
solenoid, immersed in liquid Helium was used to ap¬ 
ply magnetic fields perpendicular to the sample plane. 
The back gate characteristics of four different devices 
(see Fig. 2(b)) with different Bi 2 Se 3 crystal thickness are 
shown in Fig. 2(c) which depicts the four-terminal con¬ 
ductivity a as function of applied back gate voltage V g . 
For the 28 and 30 nm thick Bi 2 Se 3 samples almost no 
gate tunability is observed, which is in contrast to the 
two thinner (12 nm and 16 nm thick) samples where a 
can be tuned by a factor of around 2. In none of our sam- 
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Figure 2. (color online) (a) Schematic illustrations of the 
sample geometry of contacted Bi 2 Se 3 flakes. In the left panel 
we highlight the width (VF), length (L) and source-drain dis¬ 
tance ( Lsd ) of a contacted sample. The source-drain dis¬ 
tance in our samples is in the range of 14 — 19 pm. Right 
panel shows a cross-section of our samples highlighting the 
flake with thickness ( t ) resting on SiC >2 with back gate, (b) 
AFM images of the four investigated devices contacted with 
electrical contacts for Hall effect and magneto-resistance mea¬ 
surements. The scale bars are all 5 pm. For more details on 
the geometry please see Table [I] (c) Conductivity as function 
of back gate voltage for all four samples, (d) Two-dimensional 
carrier density n of the different samples, extracted from the 
Hall resistance vs back gate voltages, (e) Resistivity as func¬ 
tion of carrier density for all devices. The inset shows the 
carrier density at zero back gate voltage as function of Bi 2 Se 3 
flake thickness. 



B field [T] 



Figure 3. (color online) (a) Magneto-resistance A p = p — p° 
[where p° = p(B = 0 T)] at zero back gate voltage for the four 
devices in Fig. 2. The inset shows WAL dips at small B fields, 
(b) Ap as function of F>-field for the 16 nm thick sample for 
various back gate voltages (see labels). The inset shows the 
slopes of the linear magneto-resistance at high magnetic fields 
(see dashed lines in main panel) vs 1/n 2 . 


thickness of the SiC >2 gate oxide with a dielectric constant 
of e r « 4. From the vertical offsets of n in Fig. 2(d) we 
can estimate the average bulk charge carrier density. The 
inset in Fig. 2(e) shows the 2D carrier density n at zero 
gate voltage as function of the Bi 2 Se 3 crystal thickness. 
From the slope of this linear dependence (dotted line) we 
extract a 3D bulk carrier concentration of 1.5 xlO 19 cm -3 
in our VSS Bi 2 Se 3 crystals. Finally, we plot the resistiv¬ 
ity p = l/cr of all four devices as function of n (main 
panel of Fig. 2(e)). The observed overall trend highlights 
a consistent carrier density dependency on the measured 
resistivity of all measured devices. The increasing gate 
tunability at lower carrier densities might be either con¬ 
nected (i) to the linear density of states of the 2D surface 
states or (ii) to a reduction of the 3D bulk density of 
states (or diffusion constant) at lower Fermi energy val¬ 
ues. By assuming that surface states dominate the gate 
voltage dependence at low carrier densities we can esti¬ 
mate the carrier mobility from the nearly linear increase 
of a as function of V g (see dashed lines in Fig. 2(c)). 
From our data we extract respective surface carrier mo¬ 
bilities of p& 1600 cm 2 (Vs) -1 . As the condition pB > 1 
can be reached with experimentally accessible magnetic 
fields, Shubnikov-de Haas oscillations should become vis¬ 
ible for 5-fields larger than 5 T. We note, however, that 
we do not observe any Shubnikov-de Haas oscillations 
(see magneto-transport measurements below). We there- 


ples we observe an ambipolar transport behavior, which 
is a first indication that very high n-doping of Bi 2 Se 3 is 
present in all our devices. 

We performed Hall effect measurements to deter¬ 
mine the charge carrier densities. The extracted two- 
dimensional (2D) electron density, n, which varies lin¬ 
early with V g is shown in Fig. 2(d) for all devices. 
The slope for the three thinner samples of 7.5 — 
8.2 x 10 10 cm -2 /V (see table 1) is in reasonably good 
agreement with the geometrical gate lever arm a g = 
6o6 r /(\e\d) « 7.2 xlO 10 cm -2 /V, where d = 285 nm is the 


device # 

t (nm) 

L (pm) 

W (pm) 

Oi g (cm 2 V 1 ) 

1 

12 

3.4 

2.4 

7.5xl0 10 

2 

16 

3.5 

2.9 

8.2xl0 10 

3 

28 

4.7 

10.6 

7.5xl0 10 

4 

30 

3.5 

11.8 

9.8xl0 10 


Table I. Geometrical dimensions of the four devices discussed. 
The dimensions are defined as in Fig. [2ja) and were measured 
using an AFM. The respective gate lever arms a g are also 
included. 
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Figure 4. (color online) Weak antilocalization (WAL) peak in conductivity as function of B field (solid lines) and fits 
according the HLN model (dashed lines) for two different devices, (b), (c) Extracted fitting values for a (panel b) and l$ (panel 
c) are shown as a function of the charge carrier density for three different devices, (d) Phase coherence length as function of 
the dimensionless conductivity g D . Here, the solid line displays the theoretical result by AAK [41]. The dashed line highlights 
the result modified due to a finite electron spin-flip scattering time r s f (see text for details). 


fore conclude that is very likely that transport in our 
samples is dominated by bulk transport, where the gate 
tunability originates from a Fermi level dependent 3D 
bulk density of states or diffusion constant. 

For gaining more insights on the separation of bulk and 
surface transport we performed four-terminal magneto¬ 
resistance measurements (see Fig. 3(a)). The two most 
prominent features in our magneto-transport data are 
(i) a LMR at high magnetic fields [421445] and (ii) a rea¬ 
sonably strong WAL dip at low magnetic fields (below 
1 T) [44l 1461449] . as shown by the inset of Fig. 3(a). 

Above B = 5 — 6 T all four devices exhibit LMR. 
Interestingly, the strength of the LMR does not solely 
depend on the sample thickness nor the total charge car¬ 
rier density. We observe that the thinnest and thick¬ 
est Bi 2 Se 3 crystal exhibit similar LMR slopes (see black 
and orange curves in Fig. 3(a)), whereas the other two 
devices (see red and blue curves in Fig. 3(a)) show a 
significantly larger slope of the LMR. However, within 
a single device we find a systematic carrier density de¬ 
pendence of the slope of the LMR at large magnetic 
fields (see Fig. 3(b)). A more detailed analysis of the 
carrier-density dependent LMR slope shows that, inter¬ 
estingly, Ap/AB changes linearly as function of n ~ 2 (see 
inset in Fig. 3(b)). This dependence is in agreement 
with the generic quantum description of galvanomagnetic 
phenomena by Abrikosov 50452]. leading to p oc B/n 2 . 
However, although for the investigated B field range the 
required condition p,B > 1 might be fulfilled, we are cer¬ 
tainly not in the extreme quantum limit where only the 
lowest Landau level is filled. Moreover, it should be noted 
that the fit shown in the inset of Fig. 3(b) does not cross 
the origin. All these bring us to the conclusion that the 
LMR in our devices is rather dominated by the classi¬ 
cal linear magneto-resistance [52] due to bulk inhomo¬ 
geneities and defects in the B^Ses crystals which may 
also explain the high bulk carrier densities. 

These findings are in contrast to the WAL, which ex¬ 
hibits a clear crystal thickness dependence (see inset of 
Fig. 3(a)) and which is therefore - also in agreement with 


literature [48] 49] 53] - most likely a better fingerprint for 
surface state transport. WAL signatures are indeed in¬ 
herent to the 2D states of TIs 0123 HU As it is governed 
by quantum mechanical interference, a detailed investi¬ 
gation of resulting corrections to the conductance allows 
to learn more about phase coherent transport properties 
in these materials. Indeed, WAL in TIs has already been 
studied in great detail [54457] and it has been shown that 
the so-called Hikami-Larkin-Nagaoka (HLN) model [58] 
can be used to fit the WAL corrections at low B fields. 
Within the HLN model the conductivity correction is ex¬ 
pressed as 


A a = a ( B ) — a (0) = 


2 r 


ae 
27 x 2 h 


In 


h 


4BeL 



-Ml 


( 1 ) 


where T is the digamma function and is the phase 
coherence length. The value of the amplitude a is ex¬ 
pected to be —1/2 for perfect WAL in a two-dimensional 
system. For an ideal 3D TI with two independent and 
decoupled 2D surfaces the expected value for the total 
WAL amplitude is therefore a = — 1 [49]. For fitting 
our data, the symmetric and anti-symmetric part of the 
overall conductivity were separated. This is necessary 
considering the imperfect Hall bar geometry due to the 
etch-free sample fabrication process. Fig. 4(a) shows 
WAL data, fitted with the HLN model given by Eq. 0 
for the symmetric part of the data with an additional 
term for quadratic magneto-resistance at low magnetic 
fields /3B 2 . The values for a and as extracted from 
the fits are shown in Figs. 4(b) and 4(c). 


For the two thinner samples (orange and red data in 
Fig. 4(b)) a is gate tunable around a value of —1/2. This 
indicates that either the surface states are strongly cou¬ 
pled via the highly conductive bulk or that transport is 
purely dominated by the bulk. For the 28 nm thick sam¬ 
ple (blue data in Figs. 4(b), 4(c)) no gate dependence 
is observed indicating the dominance of a bulk conduc¬ 
tion channel with a Fermi level in a regime with constant 
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Figure 5. (color online) (a) Broadening of the WAL peak with increasing temperature for the 16 nm thick Bi 2 Se 3 sample, 
(b) Dependence of the parameter a as a function of inverse temperature for different back grate voltages. (c),(d) and (e) 
Dependence of the phase coherence length l$ as a function of 1 /\fT for three different back gate voltages. The solid and dashed 
lines resemble the same theoretical models as in Fig. 4(d). 


3D bulk density of state which suppresses any gate tun- 
ability. The thickest sample (30 nm) does not show a 
distinct WAL peak and is hence excluded from our WAL 
analysis. A similar trend is also observed for the phase 
coherence length which increases with larger sample 
thickness and increasing charge carrier density. 

Interestingly, such a gate-tunable phase-coherence 
length - also observed by other groups EEIES] - is in qual¬ 
itative agreement with a scattering mechanism based on 
electron-electron interactions as predicted by Altshuler- 
Aronov-Khmelnitsky (AAK) for a two-dimensional sys¬ 
tem J4TJ EH] : 

4 = % n (4 m*k B T lng D )~ 1/2 , (2) 

where ks is the Boltzmann constant, ra* is the effective 
mass and g u = ah/e 2 is the dimensionless conductivity, 
which can be directly extracted from the measured con¬ 
ductivity. Apart from the small logarithmic correction 
(only becoming important for very small conductivities) 
the phase coherence length is a linear function of g n . By 
plotting the experimentally extracted 1$ as function of g u 
(see Fig. 4(d)) we indeed can confirm this nearly linear 
dependence. Furthermore, by assuming a bulk carrier ef¬ 
fective mass of m* = 0.15m e m , we obtain the solid line 
in Fig. 4(d) without any further adjustable parameter. 
These values are a factor of 2 — 3 larger than the values 
of extracted from our WAL measurements, meaning 
that there must be some corrections to the effective mass 
or (more likely) additional sources for dephasing. 

This becomes even more apparent when investigating 
the temperature dependence of the WAL and the ex¬ 
tracted at different carrier densities, as shown in Fig. 5. 
In Fig. 5(a) we show the WAL peak for the 16 nm thick 
sample at different temperatures, highlighting its disap¬ 
pearing at elevated temperatures. The peak at small 
magnetic fields slowly decreases in amplitude completely 
disappears at 50 K. By fitting again the HNL model to 
our data we extract the temperature dependent a val¬ 
ues and phase-coherence lengths (Figs. 5(b)-5(e)). The 


prefactor a changes towards zero for increasing temper¬ 
ature, i.e. decreasing 1/T, as seen in Fig. 5(b). More 
insights can be gained when investigating the tempera¬ 
ture dependence of the phase coherence length. In order 
to highlight the expected temperature dependence given 
by Eq. 2 we plot 1$ as function of 1 /y/T in Figs. 5(c) to 
5(e) . 


Similar to Fig. 4(d), the solid lines illustrate the esti¬ 
mates for obtained from the AAK theory (/^ oc 1 /y/T) 
for different back gate voltages, i.e. carrier densities 
which are color coded in Figs. 5(c)-5(e). Indeed, above T 
= 7 K (below 1/y/r ~ 0.4 K -1 / 2 ), the experimentally ex¬ 
tracted values are inversely proportional to the square 
root of the temperature. However, at lower temperatures, 
shows a carrier density dependent saturation behavior, 
which can not be explained by electron-electron interac¬ 
tion limiting the phase coherence length. To account for 
these discrepancies, we follow Ref. [8U] and include an 
additional inelastic electron spin-flip scattering time, r s f. 
Thus the phase coherence time = l^/D will be lim¬ 
ited by r s f at low temperatures. This leads to an overall 
scattering rate which is the sum of spin flip and the AAK 
decoherence rate, r^ 1 = + &bT In # n / (hg n ). We use 

this expression to estimate = yjDr By assuming 
a linear carrier density dependency of the spin flip scat¬ 
tering time r s f = f3n with (3 = 1.2 x 10 -24 cm 2 s, we 
obtain good agreement with all our experimental data 
(see dashed lines in Fig. 4(d) and Figs. 5(c) to (e)). The 
extracted r s f values are on the order of 10 ps. This is 
a first experimental estimate of the spin flip scattering 
time in B^Ses from WAL data. The short timescale 
most likely results from the strong spin-orbit interaction 
in this material class [62]. We emphasize that the r s value 
is not attributed to surface but rather to bulk transport 
properties. 
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IV. CONCLUSION 

In conclusion, we used a catalyst-free vapor solid syn¬ 
thesis growth method for obtaining well-shaped single- 
crystalline Bi 2 Se 3 flakes with thicknesses in the range of 
a few nanometers. We performed low-temperature trans¬ 
port measurements on such B^Ses flakes with different 
layer thicknesses, resulting in different (two-dimensional) 
doping values. From magneto-transport measurements 
we extract information on the linear magneto-resistance 
as well as on phase coherent transport properties. In par¬ 
ticular from weak antilocalization measurements we gain 
detailed insights on the phase-coherence length in bulk 


transport. We observe that the phase-coherence length 
linearly depends on both the conductivity and electron 
density. Its values are close to the values imposed by 
electron-electron interaction but limited by spin-flip scat¬ 
tering at the lowest temperatures. 
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